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Abstract. We propose a new method of the analytical computation of the spectral 
dimension which is based on the equivalence of the random walk and the q-state Potts 
model with non-zero magnetic field in the limit (7 — > 0. Calculating the critical exponent 
of the magnetization S of this model on the dynamically triangulated random surface 
by means of a matrix model technique we obtain that the spectral dimension of this 
surface is equal to two. 
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1. Introduction 

A diffusion in the dynamic medium is an important problem due to its wealth of 
applications in many branches of physics. Examples include the diffusion through fluid 
membranes [HE], the diffusion in the presence of two-dimensional quantum gravity [31 H] 
and others. The main characteristic of the diffusion (random walk) is the probability 
that the particle returns to the starting site at time t. At large times one expects 
this probability to decay as t^*^"/^, where dg is the spectral dimension of the underlying 
geometry. In the case of Z'^ we have dg = d, which simply gives the dimension of the 
regular lattice. It is interesting how the fractal structure of random lattices affects the 
spectral dimension. Despite the fact that the spectral dimension of many random graphs 
has been calculated analytically, in such branched polymers or generic trees 

(is = 4/3 [3 E], non-generic trees or multi-critical branched polymers with k phases [7j: 



random combs ds = (4 — 6)/2, 6 < 2 [8], where 6 is a power law exponent for the length of 
the tooth of the random comb, there is no theoretical derivation of the spectral dimension 
of random lattices of a given, for instance planar or toroidal, topology. Dynamical 
triangulations, dual to random lattices, arise as a discretization of the integral over the 
metrics of some smooth two-dimensional manifold [9]. Each triangulation (see Fig. 1) is 
in one-to-one correspondence with a vacuum diagram of some N x N Hermitian matrix 
model. In the large limit only lattices with planar topology survive [TOl ITT] . 

We show that the random walk is exactly equivalent to the g-state Potts model 
with non-zero magnetic field taking in the limit q ^ [12]. On a random lattice this 
model is defined by the partition function: 



where is the adjacency matrix of the graph, the upper index n is the number of 
vertices in the graph, indices i,j = 1 . . .n enumerate vertices, spin variables cXi associated 
with the vertex i take q different values (colours) enumerating independent components 
of the spin, YIg represents the sum over all configurations of graphs in the ensemble, 
represents the sum over all configurations of spins, H represents the magnetic field 
normalized by the temperature, < a < 1 is an auxiliary parameter which is essential in 
the limit of small q, (3 is the product of the inverse temperature and the coupling constant 
of spins. The limit g ^ can be better understood through the cluster representation 
[T3] of the model (11. ip . It also can be defined in terms of the tree-like percolations 
(spanning forests) on a random graph [TT]. In this representation the partition function 
( II. ip for zero magnetic field = is given by 



2k + 2 
2k + I 



A; = 2,3,4..., 






Gn trees 
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where Yltrees trees spanning the lattice b{tree) is the number of 

bonds in a given tree on the lattice Gn and the constant B is connected to (3 in (11. II) 
through the equality = 1 + q'^B. 

The behaviour of the return probability, the conditional probability and the 
square displacement of the random walk at large times can be found by computing 
critical exponents of the magnetization 6, the two-point correlation function rj and the 
correlation length u of the spin model (11. ip correspondingly [T5] . 

The Potts model on a random lattice [HI Il7j belongs to the long list of exactly 
solvable models which could be reformulated as matrix models. The list includes the 
Ising model [IHl[19], bond-percolations [20], tree-percolations [21], the 0{n) vector model 
[221 [23], dilute Potts model [21] and many others [251 126]. The limit g of the 
Potts model which is relevant for our consideration had been solved by the saddle point 
technique [21], by the loop equation technique [IT] and, recently, by the combinatorial 
method [28]. We generalize this results to the case with non-zero magnetic field which 
breaks the symmetry of the model. In this scenario there will be two different saddle 
points. However we shall show how it can be simplified in the limit g — 0. 




Fig.l: Duality between fat graphs arising from the perturbation expansion of the 0^ of 
one-matrix model and triangulated surface. Each (p^ vertex of the fat graph corresponds 
to the triangle. Gluing together triangles edge to edge is equivalent to the Wick's 
contraction of (p matrices. 

We derive the exact result for the spectral dimension of a random surface and 
show that in the case of non-zero field there are two phases. One phase corresponds 
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to Hamiltonian walks [22] or branched polymer phase dg = 4/3 {'jstr = ~1) another is 
dilute or pure gravity phase ds = 2 {•jstr = 

This paper is organized as follows. In section 2 we establish the equivalence of the 
random walk and the Potts model. We probe this equivalence computing the critical 
exponent of the Potts model magnetization on the Bethe lattice and comparing it to the 
spectral dimension of the corresponding lattice. In section 3 we reformulate the Potts 
model on a random surface as a random multi-matrix model. In section 4 we find the 
solution of this model and construct the phase diagram. In section 5 we compare our 
results with numerical simulations and present our conclusions. 

2. Random walks and the Potts model 

A lattice is a set of n vertices connected by links. It is uniquely defined by the n x n 
adjacency matrix gIj\ whose entries are G[f = if there is no link between i and j 
and Gj^j ^ > otherwise. The coordination number or the degree of a vertex is defined 
as the total number of links connected to it. Consider the ensemble of lattices G„ with 
n vertices and the random walk jumping on the sites of the such lattices. At each time 
step the walker with equal probability must jump to the nearest-neighbour site. This 
process is independent of what can happen to the lattice bonds. We consider that lattice 
changes its configuration by choosing new one from G„ at random every time step. A 
jump can occur only if the sites are connected by the bonds at the time the walker 
attempts to jump. Let G^^\k) denote the adjacency matrix at time step k. Suppose 
that the walker starts at time /c = at the site 0. Given a particular bond history 

G = G(A:) = {G;;)(0)...G;f(fc)} (2.1) 

let Pi{k\ G) be the probability to find the random walk at the vertex i after time t and 
Pi(0) = 5o,i- The master equation for conditional probabilities of the random walk on 
the dynamical lattice of n sites can be cast into the following system of equations 

k 

p.{k + 1;G) = Y, {G'^\k)p,{k- G) + [l- Gl;\k)]p,{k; G)) , (2.2) 
i=i 

where the lower case indices i,j enumerate vertices and k counts the number of jumps 
made by the walker. The random walk is non-Markovian because jumps depend on 
bond histories G. 

It was shown numerically in [30] that conditional probabilities Pi{k; G) in the long 
time limit k +oo did not depend on the particular bond history and approached 
some average value Pi{k). This universal behaviour is very similar to the behaviour of 
the spin system on the fluctuating lattice. 

We show that there is an explicit correspondence between the random walk and 
the Potts model. Firstly, we consider the simplest case when the bond configuration 
is static. Then conditional probabilities Pi{k) become Markovian because all jumps are 
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independent of earlier events. Then the master equation (12.21) simphfies and after taking 
the continuum time hmit it can be written as the system of differential equations 

P« = E^^^feW-P^W)' (2.3) 
j 

where Gij is the adjacency matrix of some fixed lattice. After the Laplace transform 

r+oo 

P,{H) = / p,{t)e-'''dt (2.4) 
the master equation (12. 3p becomes the system of linear equations for the quantities 

Pi = Pi{Hy. 

^ij^i ~ ^o,i Lij = HSij — Gij, (2.5) 

where Lij is the Laplacian of a graph with the the adjacency matrix Gij. This system 
can be easily solved by inverting the Laplacian. The return probability is defined by 
the determinant of the Laplacian 

n on 

On the other hand, it was rigorously proven in [31] that the determinant of the Laplacian 
is the sum over all spanning forests on a lattice. The forest, including / trees which span 
mi, i = 1 . . .1 vertices correspondingly, gains the weight if' Yli '^ii where mi = n and 
n is the total number of vertices in a lattice. Thus it can be interpreted as the cluster 
representation of the partition function of g ^ Potts model with non-zero magnetic 
field. The return probability is the magnetization M and it scales as M ~ H^^^ at the 
critical point. The spectral dimension is defined by 

4 = 2 (^1 + . (2.7) 

We consider a generalization of the above result to the case of a random lattice. 
The key conjecture is 



lim 



{2.1 



where Zn \(3, H) is the partition function (II. ip . To see if it is true we consider the graph 
associated with the path of the walk on a random lattice. The edges of the graph are 
those bonds on a lattice that the random walk crossed for the first time on its path. It 
is easy to see that this graph is a tree which is called the forward tree (see Fig. 2). Thus 
the partition function of the random walk which is the sum over all possible paths of 
the random walk is dual to the partition function of all trees on a random lattice [Ij. 

Our conjecture (12.81) is actually a powerful tool for computing exponents of the 
random walk on different graphs. It allows to employ a critical phenomena technique 
to the random walk problem which is much broader than methods of a direct solution 
of the equation (12.21) . 
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Fig. 2: Forward tree is pictured for the path of the random walk on the square regular 
lattice. The walker starts at the lower left corner and ends at the upper right corner. 
Jumps through the diagonal are also allowed. Arrows denote the direction of jumps. 
Edges, belonging to the path of the walker, are shown by broken lines. The edge belongs 
to the forward tree if and only if the walker, jumping through this edge, gets to the vertex 
it had never visited before. All such edges are shown by solid lines. 

As an example we consider the Bethe lattice with coordination number z (see Fig. 
3). We put Potts spins at each site of the lattice. The solution of this model is well- 
known [32]. In the thermodynamic limit the magnetization of the central site of the 
lattice is 



M = coth.[{H - zs)/ 21 



[2.9) 



where s is a parameter defined hj x = e"^ and x is the fixed point of the recurrence 
relation which for the case of g — > can be written as: 

, = £^±e^l)£ri. (2.10) 



Recasting the exponent of the magnetic field from (12.101) we find that 



„2-l' 



2 + x 



X- ^— -. (2.11) 

e^x — I 

Now we can write H as Taylor series of the small parameter s. Up to two leading terms 
the expansion proceeds 

H = {z-2)s + {e^ -1)-^^ + .... (2.12) 

On the other hand from the expansion of the magnetization (12.91) one has M ~ s^^. 
Treating s as a function of M we obtain 

H = {z- 2)M-^ + {e^ - ly^M^^ + ... (2.13) 

Using the scaling hypothesis H = fs{M^^/^) one has 

H = M-^f,{M). (2.14) 
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It gives the values for critical exponents 6 and /3 of the g-state Potts model g — > on 
the Bethe lattice 6 = —2, jS = —1. Knowing 6 and using formula (12.71) we derive that 
ds = 1 in agreement with the result [331 El] • 





Fig. 3: Bethe lattice with coordination number z = 3. 



dMi... dMg exp{N tr[c ^ M^M, - ^ + ^^Mf + ^ ^Mf)]), 

i^i 1=1 1=2 



3. Correspondence with the Matrix Model 

From now on we restrict ourselves to the ensemble of random lattices with coordination 
number 3 and the topology of the sphere. Consider the q-matrix model defined by the 
partition function 

(3.1) 

where Mj, i = 1 . . .q are N x N Hermitian matrices, H is the magnetic field, < a < 1 
and 

c=l/(e^ + g-2), (3.2) 

where — 1 = q"'B. This model generalizes the model of [171 |2T] to the case of non-zero 
magnetic field. 

We note that all matrices in ( 13. ip are coupled to each other. Physically each matrix 
Mi represents one component of spin. It can be shown that the propagator is 

(l-c(g-2))/c,^ = j 
(c+l)(l-c(g-l)) i l,ty^J 
where (. . .)o denotes the Gaussian average {g = 0). Using Feynman diagrammatic 
expansion one would get that the free energy corresponding to (13.11) is equal to the 
generating function: 

+00 y \ n 

"^0^ V(c+l)(l-c(g-l))y 

where Zn\p,H) is the partition function (II. ip . 

After the change of variables Mj — * Mj(2(l + c))^^^"^ we have 

,2 1 9 nf^Hq"" 1 



(trM.M,)o = iV ^^ , ^^^^ ^ ^ { ^33^ 



/wZ 1 — Hq ^ — 

dM,... dM,exp{NtT[^Y' - 3 E^' + + E f ^')])' 



(3.5) 
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where g = g{2{l + c))~^/^, h'^ = c/(l + c) and Y = Mi + . . . + Mg. By introducing in 
f l3.5p new Gaussian-distributed random matrix variable X we replace the first term in 
the exponent by the matrix polynomial linear in Mf. 

dX J] dM, exp(iV tr[-XV2 + /,X ^ M, - ^ + ^^Mj^ + ^ |Mf)]). (3.6) 

i 1=1 1=1 1=2 

We want to express the integral over matrices (13. 6p by the integral over eigenvalues. 
As it was demonstrated in [35] the integral over the matrix in the external field can be 
reduced to the integral over eigenvalues by the following formula: 

[ dMexp{NtT[-My2 + MX]) = j JJ dmi^y^ exp(A^[-m,V2 + m,x,]) (3.7) 

Using (13. 7p and noticing that all integrals over Mj in (13. 6p are similar, the partition 
function can be rewritten as: 



where 



j JJrfa;iA(x)2-«exp 1x^x^2] Q+{x)Q^ 
i=i \ i=i J 



r ^ ( ^ 1 oe-^^" \ 

e+(x) = y n dmiA{m) exp f iV J][/iXimi - -m^ + ^^——mf] 1 , (3.9) 

N [ ^ I ' \ 

JJ^ (imjA(m) exp j A^^^[/ia;jmj ml H — m\] j . (3.10) 

t=i \ i=i 2 3 J 



xy-\ (3. 



e_(s) = 

In the present paper we will not give the general solution of (13.80 . For our purposes it is 
enough to find the partition function when g — 0. In the absence of the magnetic field 
H = functionals Q+{x) and 9_(a;) are equal to each other. Hence the value of this 
functionals will be governed by the same saddle point equation as shown in [201 121] . By 
noticing that h"^ = l/q°'B and after the change of variables mj mi/ y/q^-B one has 

e+(a;) = I dmAim) exp(^ ^[ximt - ^mf - ^ — m^]), (3.11) 
^ ^ i=i 

where G = g{\/2q"'B)~^ . When q is small the prefactor in the exponent becomes large 
and the steepest descent method can be used to compute (13. lip . Unlike the usual large- 
N limit the contribution from the Van-der-Monde determinant will be small and can be 
neglected. To the leading order we have 

e+(x) = exp f Y}x,u, - -u] - ^^?] j (3.12) 
Similarly we have 
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where Ui and Vi are given by the saddle point condition: 

Xi = Ui + Ge^^^M- Xi = Vi + Gvf. (3.14) 

In the hmit of small q one can express the solution Ui in terms of Vi as perturbation 
series. Choosing the ansatz 

Ui = Vi + evf^ e = Hq" (3.15) 

and expanding the exponent e^'^" we have 

By doing this the ratio of (13.121) and (13.131) is significantly simplified: 



Q+(x) I NGH ^ 



exp 



E-n (3-17) 



e_(x) " \ w 

It follows that the partition function is 

Z = j n^x.A(x)^exp {^^jl^V^ (-^E-.') (3-18) 
After changing variables from x to w and using (I3.14p the integral (13.181) becomes 

Z = j dvi\{{l + G{vi + Vj))/\{vfeM-NY,y{'"i)l (3-19) 

n^)=Q(^ + Gx^f + ^). (3.20) 

Again, we apply the steepest descent method but now with respect to large A^. The 
saddle point equation is 

2 E + E N = ^^'(^0- (3.21) 

The distribution of the eigenvalues becomes continuous with density p{x) = 
(l/N) ^5(x — Vi). We restrict ourselves to the one-cut case where all eigenvalues Vi 
belong to the support consisting of one interval [a,b], ah > 0. The equation for the 
density is the integral equation: 

P{y)dy ^ I p{y)dy ^ l^,^^^^ ^3 ^2) 



X — y J 1/G + x + y 2 
where P denotes the principal value of the integral. Introducing the trace of the resolvent 

/ N 1 1 [ p{y)dy 

co'(x) = — tr— = / (3.23) 

^ ^ A^ M-x J y-x ^ ' 

(I3.22P can be equivalently written as 

w(x + iO) + uj{x - iO) + 2lj(-x - l/G) = ]^l\x) (3.24) 
The integral equation that governs the eigenvalue density is 

+ (3.25) 



[x - y){\ ^ G[x ^ y)) 2' ' 2B\^2Gx 
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The equation fl3.25p represents the Riemann-Hilbert problem. The solution can be found 
ESI: 



^2w,2_ 2\ii/2 fiy)dy 



p{x - 1/2G) = -[{x' - a'){h' - x')Y'' / — , iT, (4.1) 

where 

f{x) = Gx^ + HGx^/{2B) - (1/4 + H/{2B))x + H/{8BG), (4.2) 
supplied with additional condition: 

/' im = (4 3) 

and with the normalization condition: 

-b 

dxp{x) = 1. (4.4) 



The ends of the interval [a, b] are functions of G, B, H, to be determined from 
transcendental equations (14. 3 p and (14.40 on a = a{G, B, H) and b = b{G, B, H). 
We will need the first equation which can be resolved as 



^TT / 1 „A HGa f b'^\ fl H\7i H , u \ , , 



where K{x), E{x) are elliptic integrals of the first and second kind respectively. We used 
the following integrals: 



liy'- 




^2)]l/2 




ydy 




liy'- 




^2)]l/2 




y^dy 




[{y'- 


a2)(62- 


^2)]l/2 




y^dy 





a2 



f, (4.7) 



aEi^l--], (4.8) 

TT / 1 



Then the critical behaviour of the partition function (13. ip is obtained by taking the 
double scaling limit B Be (infinite random surface) and G ^ Gc (infinite trees). The 
latter occurs near the upper edge b of the support. We have be = —\/2Gc- It can be 
recast in the following form 

w = l±l/h,h>l l = cos(7r//i) h = 2,A... (4.10) 

Then the critical density is 

p{x) ^{g- g.)(b - xf-^"' + {b- xf^^"' (4.11) 
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and string susceptibility is ■jstr = ~l/h. The solution is singular over [a—l/2G, 6— 1/2G] 
and [—b — l/2G,—a — 1/2G]. Branching points coincide when a = [21j. It can be 
shown that the Potts model is in the critical point if this condition is satisfied. Thus 
there are two phases of the model: dilute and dense. In the dilute phase the problem is 
equivalent to the (-2)-dimensional dynamical triangulated surface. The critical exponent 
of the string susceptibility is •jstr = — 1- In the dense phase defined by the condition 
b = —1/2G one has •jstr = "~l/2. From the equation (14.51) we can see that 

db 1 , , 

Finally, let us compute the magnetization M = —dF/dH. The free energy is 

cb nb nb 

dxp{x)V{x)+ I I dxdyp{x) p{y)^n\x — y\ + + G{x + y))]. (4-13) 



It is linear in H for the small magnetic field. Taking the derivative one can obtain that 
3b 

M ~ TT^iGiH + G^H^) = [Gi + G2H) = fs{H). (4.14) 
It means that 1/5 = 0. Substituting the exponent 6 in (12. 7p one has 

4 = 2. (4.15) 



5. Conclusion 



In conclusion, we have demonstrated an alternative derivation of the spectral dimension 
of the random surface. We note that our results match the numerical simulations [3^. It 
should also coincide with the KPZ [38] result for the conformal field theory with central 
charge c = —2 coupled to the two-dimensional gravity. The result can be generalized 
to topologies with higher genus by the DDK [39l HQ] formula. We notice that it should 
not change the value of the spectral dimension ds- The only sensitive exponent to DDK 
is the string susceptibility 7str- 
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